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ABSTRACT 

Conventional treatment of cold dark matter halos employs the Navarro- Frenk- 
White (NFW) profile with a maximum radius set at r = r 2 oo, where the enclosed 
matter has an overdensity of 200 times the critical density. The choice of r = r 20 o 
is somewhat arbitrary. It is not the collapsed (virial) radius, but does give r ~ 1 
Mpc for rich clusters, which is a typical X-ray size. Weak lensing measurements, 
however, reveal halo radii well in excess of r 2 oo- Is there a surface that places an 
absolute limit on the extension of a halo? To answer the question, we derived 
analytically the solution for circular orbits around a mass concentration in an ex- 
panding flat Universe, to show that an outermost orbit exists at v/r — H, where 
v is the orbital speed and H is the Hubble constant. The solution, parametrized 
as r 2 , is independent of model assumptions on structure formation, and is the 
radius at which the furthest particle can be regarded as part of the bound system. 
We present observational evidence in support of dark matter halos reaching at 
least as far out as r = r 2 . An interesting consequence that emerges concerns 
the behavior of rotation curves. Near r = r 2 velocities will be biased low. As 
a result, the mass of many galaxy groups may have been underestimated. At 
r = r 2 there is an abrupt cutoff in the curve, irrespective of the halo profile. An 
important cosmological test can therefore be performed if velocity disperion data 
are available out to 10 Mpc radii for nearby clusters (less at higher redshifts). 
For Virgo it appears that there is no such cutoff. 

1. Introduction 

In the theory of large scale structure formation, a pillar of the standard cosmological 
model, one important question that has not been answered to satisfaction concerns where 
exactly, viz., at what threshold of overdensity, is the boundary of a mass concentration (or 
clump) at a given redshift. The question is pertinent to any effort in obtaining a reliable 
theoretical estimate of the total mass and spatial extent of the 'dark halos' in galaxies 
and clusters. The reason has to do with the availability, from numerical hydrodynamic 



- 2 - 



simulations (Navarro, Frenk, and White 1995, 1996, 1997; but see also the even earlier work 
of Dubinski & Carlberg 1991), of a 'universal' profile (known as NFW profile) for the clumps, 
of the form 

Pm(r) = ^p— 1 

for the matter density distribution, where in Eq. (1) f denotes a physical radius (i.e. an 
invariant for bound structures), r s is a constant scale radius, 5 C is an overdensity param- 
eter, and p c is the critical density. Since at the outer radii r ^> r s the density scales as 
Pmij') ~ 1/r 3 , and there is no further change of functional form with increasing distance, 
i.e. hierarchical structure formation codes do not seem to reveal the surface radius of a 
clump, the total integrated mass is divergent unless an upper radius limit, (or cutoff) R , is 
'manually' assigned. 

Conventionally this limit is set at R = r 20 o, sometimes referred also to as the virial 
radius, defined as the radius at which the enclosed matter density is 200 times above the 
critical density 

The choice of the virial radius R is probably just a matter of convenience, e.g. for a rich 
cluster like Coma one may envisage a virial mass ~ 2 x 10 15 M Q , in which case a factor of 
200 overdensity would lead, via the equation 

800tt „ 

M 20 o = — g— Pcr 200 - (3) 

to it! pa 2.5 Mpc, which is not far from the value believed to be 'reasonable'. (Lokas & 
Mamon 2003). Historically, an analytical formula for the virial radius was derived by directly 
appealing to the virial theorem (Lahav et al 1991). For an Q m = 0.27 and Q\ = 0.73 
cosmology (e.g. Spergel et al 2006) this radius (sometimes also called the 'collapsed' radius) 
is less than half the 'turnaround' radius of Peebles (1984), and corresponds to an overdensity 
of pa 90(l+2;) 3 : see Fig. 1 of Eke, Cole, and Frenk (1996). As we shall see, the actual evidence 
points to a continuation of the NFW profile to radii well beyond this value, suggesting the 
already well accepted fact that particles need not be completely virialized before they become 
part of the clumped (or collapsed) system. 

If the virial radius lacks observational significance are there other criteria available that 
may better connect with reality? The fundamental point here is that for a pure Newtonian 
(static and infinitely old) Universe there is in principle no end to the zone of gravitational 
influence of a clump. Clearly the same is no longer true for an expanding Universe. Can an 
outermost surface for a bound structure be drawn, and be subject to scrutiny? Attempts 
to define such a surface based upon consideration of radial motion, have been made (e.g. 
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Sandage 1986). The question which cannot be answered by simply calculating radial veloci- 
ties, however, is whether bound orbits can still exist at these large distances. As it turns out, 
there is a radius within which circular orbits are sustainable, The derivation of this radius 
does not depend on any of the assumptions made about how the clump is formed. 

2. Criterion for the existence of bound orbits in an expanding Universe 

If in a flat and unperturbed FRW space-time with dimensionless expansion parameter 
a(t) and c = 1, i.e. 

ds 2 = dt 2 - a 2 (t) {dr 2 + r W) (4) 

we envisage particles located at constant physical distances from each other, or geometrical 
shapes that do not expand with the Universe, it would be more convenient to use f = ar 
rather than r as radial coordinate, because our objects of interest are tied to the ar-grid and 
not the r-grid. Moreover, if we also perform the transformation t — > t where t — t — Hf 2 /2, 
Eq. (4) will become Minkowski in form, 

ds 2 = dt 2 -(l + H 2 r 2 )(dr 2 - r 2 dVL 2 ), (5) 

apart from one correction term which is only second order in H (we assumed zero acceleration 
of the expansion by ignoring a qH 2 term; as will be evident from the complete treatment 
below, the effect of q is negligible for our purpose). 

Eq. (5) is the reason why no experiments performed on earth, or within the Milky Way, 
can directly probe the Hubble expansion. An interesting point emerges nonetheless when 
one computes the radial speed of light 

df 1 9 9 

= 1 - -H 2 f 2 , 
dt 2 

and find that to order H 2 it decreases with increasing f, symptomatic of an effective repulsive 
potential that eventually prevails at large radii where expansion might overcome gravity. At 
this stage, however, the idea is only suggestive, because we have yet to formally include the 
effect of gravity. 

If the large scale character of space-time is that of a flat FRW metric, but locally there 
is perturbation by a weak and centrally symmetric gravitational field, the line element will 
be modified to the form 

ds 2 = (1 + 2§)dt 2 - (1 - 2$)a 2 (t)(dr 2 + rW), (6) 

where <£> is the gravitational potential (McVittie 1933, see also Futumase & Sasaki 1989). 
This time the transformation to coordinates (t, f) no longer offers a great deal of simplifica- 
tion. Nevertheless, McVittie (1933) showed that circular orbits maintain their constant radii 
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r = a(t)r, so there is still some mathematical advantage to be gained in preferring r to r as 
the radial coordinate for our problem, i.e. in a (t, f) system Eq. (7) reads as 

ds 2 = [1 + 2<5>{?)]dt 2 - [1 - 2$(f)][(df - tffdt) 2 + fW]. (7) 

Although in Eq. (8) one encounters a cross term of the form Hfdfdt which is first order 
in the Hubble constant, we shall soon find out when we calculate orbital behavior that, in 
agreement with our earlier analysis, no first order effects exist. 

To compute orbits from Eq. (8), we replace each dx^ in the generic form ds 2 = g PhV dx >1 dx v 
by dx^/dr to obtain the Lagrangian for motion along the 9 = plane as 



dt x 2 



2L = [l + 2$(f)] ( — J -[l-2$(f) 



\dr dr ) ~^ \dr 



(8) 



where r is the proper time, and a factor of two was introduced on the left side by convention 
(so that the canonical momenta dL/dt = t and dL/dr = a 2 r, rather than twice these 
quantities). The three Euler-Lagrange (or geodesic) equations corresponding to t,r,9 are, 
respectively: 

i(l + 2$)i + (1 - 2$)Hf(f - Hri)} = (1 - 2$)^f£(f - Hri), (9) 
dr dt 

-i(l - 2$)(f - Hri)] = ^[i 2 + (f - #fi) 2 ] + (1 - 2$)[ffi(f - fffi) - f0 2 ], (10) 
ar dr 

2^1=0, (11) 

where the 'dot derivatives' are w.r.t. the proper time; e.g. t = dt/dr. Since Eq. (9) is a 
homogeneous function of the time derivatives, there is an immediate first integral: L is a 
constant. And to normalize r to be the proper time, we take 2L — 1. 

We are interested in the effect of expansion on circular orbits around a mass concentra- 
tion, so let us look for solutions in which (at least on time scales short compared with the 
Hubble time) f is a constant. Moreover, 6 and t are constants, the ratio between them gives 
the angular velocity: 

d9 9 

= (12) 

Note that on the right side of Eq. (10) the quantity dH/dt is either multiplied by f, which 
we have assumed is zero, or by H. Since dH/dt = —(1 + q)H 2 where q = —ad/a 2 depicts 
the (negative) acceleration of the expansion, the product HdH/ dt ~ H 3 and can be ignored 
if we are only concerned with terms of order H 2 or lower. 
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Under our assumptions, the left-hand sides of Eqs. (10) to (12) all vanish, as do the 
right-hand sides of Eqs. (10) and (12). So we are left with one equation, 



or, equivalently, 



= ^(1 - H 2 r 2 )i 2 - (1 - 2§)(H 2 ri 2 + f6 2 ), (13) 
dr 



(1 - 2$)(tf 2 + u 2 )f = (1 - H 2 ? 2 )^ = (1 - H 2 f 2 )^f (14) 

dr r 2 



where in the rightmost expression M denotes the excess mass within radius f, after the 
mass contribution to that region from the mean density of matter and 'dark energy' in the 
Universe is subtracted. Solving for the angular velocity of the circular orbit, we obtain 



1 - H 2 r 2 GM 

I I 2GM ~ 

f 

Provided that GM jr <C 1 and H 2 f 2 <C 1, Eq. (16) may be recast as simply 

GM 

f 



" 2 = -r-*jM^r-H 2 - (15) 



u: 2 = — -H 2 , (16) 



i.e. the net effect of the Hubble expansion is to reduce the angular velocity w.r.t. its value 
under the scenario of pure Newtonian gravity (or more precisely a Newtonianly perturbed 
Minkowski line element). A remarkable feature of Eq. (17) is that it predicts the existence 
of an outermost bound orbit, which has a radius R given by 

G J^o^ R =H. (17 ) 

Eq. (18) yielded a radius smaller than that of Eq. (4). 

It is also possible to define the outermost orbit in terms of an overdensity criterion. If 
in Eq. (18) we write M = 4nR 3 p c i mnp /3, where p c ium P is the mean overdensity within radius 
R, the criterion will read 

3H 2 

Pclump - (18) 

For mass clumps at z = we have H = H , and the rightside of Eq. (19) becomes twice 
the critical density, 2p c . Thus we shall henceforth refer to the boundary radius of a self- 
gravitating mass clump as 

fGMY , x 

r ^{-w) • (19) 

even if the terminology is obviously loose, in the sense that for z > clumps we have H > H 
and the rightside is no longer strictly equal to 2p c . 
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3. Observational consequences: rotation curves and the total matter budget 

Although the impression one gets from a superficial perusal of the X-ray images of 
clusters of galaxies is that clusters extend to radii of 1 - 2 Mpc, commensurate with the 
value of r 20 o as defined in section 1, there are plenty of evidence pointing to the existence of 
massive extended halos in clusters and galaxies, i.e. the cutoff radius is more consistent with 
r 2 of Eq. (20) than with r 20 o- We provide examples in each case to substantiate our claim. 

3.1 Clusters 

A typical number for the lower mass range of rich clusters is M = 10 15 M Q . By Eq. 
(20) a z = cluster has boundary radius at r 2 ~ 10 Mpc, whereas the same for a z = 0.4 
cluster is r 2 ps 8.5 Mpc (we adopted h = 0.7 throughout this subsection). Thus, whether one 
'weighs' a cluster by measuring the velocity dispersion of member galaxies, or assembling 
a weak lensing shear map from background sources, the easier task is to target at higher 
redshifts: not only is a cluster's angular size smaller because of the larger distance, but 
also its limiting physical radius r 2 is reduced. Weak lensing observations did reveal a very 
extended halo for the z = 0.395 cluster CL0024+1654 (Kneib et al 2003). The cluster mass 
as quoted in this paper is M 20 o ~ 6 x 10 14 M , so that one expects r 2 > 7.2 Mpc (most 
likely much larger). Indeed, the data indicated a density profile p m (r) ~ l/r 2A continuing 
through the instrumental sensitivity limit at f ~ 5 Mpc without any sign of cutoff. 

An important test of the standard model can be conducted if rotation (velocity disper- 
sion) curves of clusters are measured out to f = r 2 . This means, by Eq. (20), that one needs 
data out to 10 Mpc radii for nearby clusters, smaller for higher redshift clusters. Assuming 
the halo profile at the outskirts has the asymptotic NFW form, then, by Eq. (17), the 
observed circular velocity will scale with radius as 

•J. = G(M ° + f M ' hrf) - Hf>, (20) 

where M c = Aiir z s p c 5 c and Mq = M 2 oo — M c lnr 2 oo- Thus, while pure Newtonian gravity (the 
first term on the right side of the equation) predicts an essentially v b s ~ 1/ Vf decline for 
r s <C f <C r 2 , for an expanding Universe as f — > r 2 there is a sharp dive towards t> b s = 0. 
If the halo profile is an isothermal sphere, where p(f) ~ 1/r 2 , we will have a flat rotation 
curve (i.e. v Q \, s = vo, a constant) when H = 0. With expansion, however, t>^ bs = v q — H 2 r 2 , 
again an abrupt cutoff at f = r 2 . 

For the Virgo cluster, the data for this test are either available or (with databases like 
the SDSS) imminently so. The fact that there is controversy over how the Local Group 
relates to the cluster implies galaxies as far out as 12-13 Mpc radii are still members of 
Virgo. Note that at these radii H r — > 1,000 km s -1 . i.e. the true velocities (due to gravity 
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alone) have to be well in excess of H r in order to maintain the rotation curve against rapid 
decline. Thus, either Virgo has a total mass ^> 10 15 M Q , or the predicted outermost orbit 
does not correspond to reality. A closer look into this problem is definitely priority task. 

3.2 Galaxies 

Weak lensing of background sources by foreground galaxies was investigated by Hoek- 
stra, Yee, & Gladders (2004), and Hoekstra et al (2005), who found clear signals out to at 
least f = r 2 oo- Based upon the lensing data alone, no definitive statements could be made 
about what lies beyond, apart from the fact that there was no indication whatsoever of r 20 o 
as representating any real cutoff radius. It is possible, nonetheless, to derive an average mat- 
ter density p g for galaxies from the Hoekstra et al observations, and to investigate whether 
the inclusion of extra matter between r 200 and r 2 would lead to a revised value for p g that is 
closer to expectation. 

According to Hoekstra, Yee, & Gladders (2004), the mass within f = r 2 oo = 139/iT 1 
kpc is M 20 o = 8.4 x 10 n /i -1 M when the data for a representative galaxy were modelled 
with the NFW profile. From Eq. (20), we see that even without any outlying matter, the 
r 2 radius for M = M 2 oo is r 2 m 1/tT 1 Mpc. Since, in this paper the scale radius of the NFW 
profile has the fitted value of r s rs 16/i" 1 kpc, we have f ^> r s for all radii r 20 o < f < r 2 , 
so that the NFW profile reduces to p s (f) = 5 c p c rf/f 3 , which can readily be integrated from 
f = r 20 o to f = r 2 . The outcome is, of course, the mass of the remaining galactic matter 
halo, viz. 



By using the best-fit 5 C parameter of 5 C = 2.4 x 10 4 as quoted in the paper, we then deduce 
that M(r 2 oo < f < r 2 ) equals ~ 85 % of Mif < r 20 o). It is clear also that in reality the halo 
boundary has radius > lh~ l Mpc because of the extra mass from the r 2 > f > r 2 oo region 
which we did not take into account when calculating r 2 . Thus we arrive at the comparison 



for NFW galaxy halo profiles, i.e. the inclusion of outlying matter would usually lead to a 
doubling of the total mass. 

In order to derive p g , it is necessary to combine galaxy mass and luminosity mea- 
surements, and to obtain a luminosity density for the same sample. The former was done in 
Hoekstra et al (2005), which still reported a representative value of M 20 o close to that quoted 
in the previous paragraph, and which also determined an average B band mass-to-light ratio 
of M 200 /L m 60 h M /L (see the bottom left plot of Figure 8 of Hoekstra et al 2005). The 
latter is to be extracted from the CNOC2 survey of Lin et al (1999), which targeted exactly 




(21) 



M(f < r 200 ) « M(r 200 < f < r 2 ) 



(22) 
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the same range of galaxies as those of Hoekstra et al (2005) and Hoekstra, Yee, & Gladders 

(2004) . From Table 3 of Lin et al (1999), the total B band luminosity density for 0.25 < z < 
0.4 galaxies (same redshift interval as that of the Hoekstra et al 2005 lens sample) is pf = 
1.2/ix 10 20 W Hz" 1 Mpc~ 3 . Since the B band solar luminosity is L% = 2.19 x 10 11 W Hz" 1 , 
we may now couple this with the aforementioned mass-to-light ratio to arrive at an average 
mass density of 

p g « 3.29 x 10 w h 2 M Mpc" 3 . (23) 

When comparing with the critical density of 6.11 x 10 n h 2 M Q Mpc~ 3 at z — 0.3, the mean 
redshift of the Hoekstra et al (2005) lenses, we see that the galaxies account for 5.4 % of the 
matter at z — 0.3. 

Is this percentage reasonable? From the matter budget analysis of Fukugita (2004) and 
Fukugita, Hogan, & Peebles (1998) emerges the picture that ~ 50 % of the matter in the 
near Universe is still 'missing', and may well reside in galaxies and their extended halos. If 
Q, m 0.27, then the expectation is f2 g ^ 13.5 %, i.e. ~ twice the percentage value as our 
5.4 %, which of course was inferred from the galaxy mass-to-light ratio where mass refers to 
M 2 oo- Given however, that we demonstrated a doubling of a galaxy's mass when the matter 
between f = r 2 oo and f = r 2 is included (i.e. Eq. (22)), the conclusion of a very extended 
halo component of baryons and dark matter fulfilling the anticipation of Fukugita (2004) 
would appear in order. 

3.3 Groups 

Groups of galaxies are the 'dark horse', in the sense that they pose a major systematic 
uncertainty to the matter budget of the near Universe, also in the direction of raising the 
fraction of matter in halos. The difference from clusters and galaxies is that, while the former 
harbors negligible fraction of Vt rn and the latter ~ 50 % of VL m if halos are included, groups 
can potentially be the refuge for a great deal more matter than either of them. Lieu & Mittaz 

(2005) analyzed the ESO survey database of 1,168 nearby groups (Ramella et al 2002) and 
found, at h — 0.7, a number density of 1.56 x 10~ 4 Mpc~ 3 and a mean virial mass of M ~ 
1.15 x 10 14 M per group. This leads, without inclusion of any extra mass that may be in 
halo extensions, to a group matter density of Vl m /2. Moreover, there is another effect. In the 
ESO survey, the peak velocity dispersion is a ~ 70 km s _1 , with a corresponding virial mass 
ofMK 10 12 - 5 M (see Figure 1 of Lieu & Mittaz 2005). The radius at which a applies is 
therefore at f GMj (2cr 2 ) = 2.66 Mpc, so that H f for h = 0.7 is not far below the circular 
velocity v = V2&, i.e. we are again in a regime where, by Eq. (17), the observed velocity 
(hence inferred mass) for these groups is biased low by the Hubble effect on the orbits. Thus, 
however one looks at them, groups have a tendency to become 'heavier', which is why some 
authors (e.g. Guimaraes, Myers, & Shanks 2005) believe that the matter content of groups 
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is a number to be reckoned with. 

In conclusion, the chief emphasis of this paper concerns the derivation of r 2 , and the 
observational evidence for its significance, including the estimate that galaxies can account 
for the missing 50 % of the WMAP matter density at low redshifts if the halo mass between 
r = r 2 oo and r = r 2 is counted. An important test of the standard model is to search for 
the decline of cluster velocity dispersions near f — r 2 . The model will be challenged if this 
decline is not seen out to radii r > r 2 , as seems to be the case for Virgo. 

We thank T.W.B. Kibble, N. Gnedin, S. DeDeo, and M. Chodorowski for helpful dis- 
cussions, especially for pointing out the absence of first order effects. 
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